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Abstract 



The manifestation of the spin-wave interaction in the low-temperature series of 
the partition function has been investigated extensively over more than seven decades 
in the case of the three-dimensional ferromagnet. Surprisingly, the same problem 
regarding ferromagnets in two spatial dimensions, to the best of our knowledge, has 
never been addressed in a systematic way so far. In the present paper the low- 
temperature properties of two-dimensional ideal ferromagnets are analyzed within the 
mo del- independent method of effective Lagrangians. The low-temperature expansion 
of the partition function is evaluated up to two-loop order and the general structure 
of this series is discussed, including the effect of a weak external magnetic field. Our 
results apply to two-dimensional ideal ferromagnets which exhibit a spontaneously 
broken spin rotation symmetry 0(3) — > 0(2) and are defined on a square, honeycomb, 
triangular or Kagome lattice. Remarkably, the spin-wave interaction only sets in at 
three-loop order. In particular, there is no interaction term of order T 3 in the low- 
temperature series for the free energy density. This is the analog of the statement 
that, in the case of three-dimensional ferromagnets, there is no interaction term of 
order T 4 in the free energy density. We also provide a careful discussion of the 
implications of the Mermin- Wagner theorem in the present context and thereby put 
our low-temperature expansions on safe grounds. 
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1 Introduction 



The question of how the spin- wave interaction in a three-dimensional ideal ferromagnet 
manifests itself in the low-temperature expansion of the partition function has a very 
long history. Dyson rigorously answered this question [lj], pointing out errors in some 
unsuccessful earlier attempts 0-[5j. After his monumental work, many researchers 
focused on how to derive Dyson's result - which is based on a fairly complicated 
mathematical formalism - with alternative methods in a simpler way Out of the 
numerous articles we would like to mention the reference by Zittarz [6], which solves 
the problem in a simple and elegant manner, as Dyson put it [3] . More recently, within 
the effective Lagrangian framework, Dyson's low-temperature series was rederived jsj 
and extended to higher orders [9J. In particular, the general structure of the low- 
temperature series of the partition function for a three-dimensional ideal ferromagnet 
was discussed in the latter reference in full detail. 

Our motivation to write the present article is based on the fact that, apart from 
some scarce papers distributed over the years, no such systematic investigation ap- 
pears to exist in the case of two-dimensional ideal ferromagnets. Above all, to the 
best of our knowledge, the effect of the spin-wave interaction on the low-temperature 
series for the partition function of two-dimensional ferromagnets has never been ad- 
dressed so far. The few available papers, all dealing with noninteracting spin waves 



10Ml9|. imply that the free energy density - for a square lattice and in the absence 



of a magnetic field - obeys the following series 

z = -fj T 2 -fj 1 T 3 + O(T 4 ). (1.1) 

However, it has never been discussed whether the spin-wave interaction already shows 
up at order T 3 or rather beyond. In other words, it remains unclear whether the above 
series referring to the ideal magnon gas is indeed complete up to order T 3 . Moreover, 
so far it has never been discussed in a systematic manner how a weak external magnetic 
field manifest itself in the above low-temperature series or how the series looks like 
for underlying geometries other than a square lattice. 

In the present work, using the model-independent and universal method of effective 
Lagrangians, we systematically evaluate the partition function of the two-dimensional 
ideal ferromagnet without restoring to any approximations. We fully take into account 
lattice anisotropics which will start manifesting themselves at order T 3 in the above 
series and thereby extend the existing results which strictly apply to the square lattice 
by considering also the honeycomb, the triangular and the Kagome lattice. We then 
show that, up to the order considered in the above series for the free energy density, 
we are dealing with noninteracting spin waves - the interaction sets in only at order 
T 4 In T. 

Even in the simple case of noninteracting spin waves, the range of validity of 
the above low-temperature series derived within the framework of modified spin-wave 
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theory 13Ml6l| - which resorts to an ad hoc assumption - has never been critically 



examined. In fact, in Ref. 20] it is stated that to systematically calculate the ther- 
modynamic properties of a two-dimensional quantum ferromagnet at low temperatures 
remains an unsolved problem of the spin-wave theory. To the best of our knowledge, 
a rigorous justification of the validity of the results obtained within the framework 
of modified spin-wave theory is indeed still lacking. In the present paper not only 
will we derive the low-temperature properties of two-dimensional ideal ferromagnets 
in a systematic manner by using effective Lagrangians, but also we will put our low- 
temperature series on a firm basis by discussing in detail the implications of the 
Mermin- Wagner theorem. 

It will also prove to be very instructive to compare the present results for the two- 
dimensional ferromagnet with those for the three-dimensional ferromagnet, adopting 
thereby a unified perspective based on symmetry. In particular, we will discuss the 
general structure of the low-temperature series for the free energy density, pointing 
out how the spin-wave interaction manifests itself in either case. 

The rest of the paper is organized as follows. In Sec. [2] we provide a brief outline of 
the effective Lagrangian method - much more detailed accounts can be found in the 



pedagogic references [2lH25j . The evaluation of the partition function up to two- loop 
order in the low-temperature expansion is presented in Sec. [3] and the thermodynamic 
properties of two-dimensional ideal ferromagnets are discussed. The range of validity 
of the low-temperature series obtained is critically examined in Sec. HI Differences 
and analogies between two-dimensional and three-dimensional ideal ferromagnets are 
discussed in Sec. [5j Finally, Sec. [6] contains our conclusions. 

We would like to mention that the systematic and model-independent effective 
Lagrangian method has successfully been applied to other condensed matter problems. 
These include antiferromagnets and ferromagnets in two and three spatial dimensions 



26M38j , as well as two-dimensional antiferromagnets which are the precursors of high- 



temperature superconductors [391-147 



Also, the correctness of the effective Lagrangian approach was demonstrated ex- 
plicitly in a recent article on an analytically solvable microscopic model for a hole- 
doped ferromagnet in 1+1 dimensions j48[, by comparing the effective theory pre- 
dictions with the microscopic calculation. Furthermore, in a series of high-accuracy 
numerical investigations of the antiferromagnetic spi n- jr quantum Heisenberg model 



on a square lattice using the loop-cluster algorithm [491-15 2|. the Monte Carlo data 
were confronted with the analytic predictions of the effective Lagrangian theory and 
the low-energy constants were extracted with permille accuracy. All these different 
tests unambiguously confirm that the effective Lagrangian technique provides a rigor- 
ous and systematic framework to investigate condensed matter systems which exhibit 
a spontaneously broken continuous symmetry. 
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2 Effective Lagrangians at Finite Temperature 



The thermodynamic properties of two-dimensional ferromagnets at low temperatures 
have been investigated before with microscopic methods, such as modified spin-wave 
theory [l3| or Schwinger Boson mean field theory 17]. The corresponding low- 



temperature series for the free energy density amounts to a power expansion in the 
parameter T / J, where J > is the exchange integral of the ferromagnetic Heisenberg 
model 

H = — J S m ■ S n — [i S n ■ H , J = const. , (2.1) 

n.n. n 

augmented by the Zeeman term which includes the magnetic field H = (0,0, H). 
The above Hamiltonian, defined on a two-dimensional lattice with purely isotropic 
exchange coupling between nearest neighbors, represents what we call ideal ferromag- 
net. 

In the present article, we will pursue quite a different approach, based on a rigorous 
symmetry analysis, which will allow us to derive the low-temperature properties of 
two-dimensional ferromagnets. One of the virtues of the effective Lagrangian method 
is that it is completely systematic and model- independent. In order not to be repeti- 
tive, here we only provide a brief introduction to the effective Lagrangian method and 
its extension to finite temperature - a rather detailed account on finite-temperature 
effective Lagrangians can be found in appendix A of Ref. [sj and in the various refer- 
ences given therein. 

Whenever a global continuous symmetry breaks down spontaneously, Goldstone 
bosons emerge as the relevant degrees of freedom at low energies. The effective La- 
grangian formulates the dynamics of the system in terms of these Goldstone bosons 



53H55I]. In the present case of the two-dimensional ferromagnet, the 0(3) spin rota- 
tion symmetry of the Heisenberg model is spontaneously broken by the ground state 
of the ferromagnet which is invariant only under the group 0(2). According to the 



nonrelativistic Goldstone theorem |56M58l|. two real magnon fields - or one physical 



magnon particle - then occur in the low-energy spectrum of the ferromagnet. 

Having identified the basic degrees of freedom of the effective theory, one system- 
atically constructs the terms appearing in the effective Lagrangian, order by order in 
a derivative expansion. The idea is rather simple: One writes down in a systematic 
way all terms which are invariant under the symmetries that have been identified in 
the underlying theory. In the present case of the Heisenberg model these symmetries 
include the spontaneously broken spin rotation symmetry 0(3), parity, time reversal, 
as well as the discrete symmetries of the square, honeycomb, triangular or Kagome 
lattice. 

The various pieces in the effective Lagrangian can be organized according to the 
number of space and time derivatives which act on the Goldstone boson fields. This is 
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what is meant by systematic: At low energies, terms in the effective Lagrangian which 
contain only a few derivatives are the dominant ones, while terms with more deriva- 
tives are suppressed. The leading-order effective Lagrangian for the ideal ferromagnet 



is of order p and takes the form 29 



C 'iff = ^ifuf + ^ HU " - \F%Wd r W . (2.2) 



The two real components of the magnon field, U a (a = 1,2) are the first two com- 
ponents of the three-dimensional unit vector U l = (U a ,U 3 ). While the derivative 
structure of the above terms is unambiguously determined by the symmetries of the 
underlying theory, the two a priori unknown low-energy constants - the spontaneous 
magnetization at zero temperature S, and the constant F - have to be determined 
by experiment, numerical simulation or comparison with the microscopic theory. 

The above Lagrangian leads to a quadratic dispersion relation 

u(k) = !k 2 + 0(\k\ 4 ) , 7 = ^, (2.3) 

obeyed by ferromagnetic magnons. This relation dictates how we have to count time 
and space derivatives in the systematic effective expansion: One time derivative (do) 
is on the same footing as two space derivatives (d r d r ), i.e., two powers of momentum 
count as only one power of energy or temperature: k 2 oc u, T. Note that, at this order, 
lattice anisotropies do not yet manifest themselves - the leading order Lagrangian 
(12.21) is invariant under continuous space rotations, although the underlying square, 
honeycomb, triangular or Kagome lattices are only invariant under discrete space 
rotations. 



As derived in [8|, the next-to- leading order Lagrangian for the ideal ferromagnet 
is of order p 4 and amounts to 

C* ff = k&irdrU*) 2 + l^drWdsU 1 ) 2 + IsAlPAU* , (2.4) 

where A denotes the Laplace operator in two spatial dimensions. The next-to-leading 
order effective Lagrangian involves the three low-energy coupling constants li,l 2 and 



The evaluation of the partition function in Refs. [8|, [9J was based on the assumption 
that the 0(3) space rotation symmetry, which is an accidental symmetry of the leading 
order effective Lagrangian, persists at higher orders in the derivative expansion. Here, 
we drop this idealization and hence also consider terms in £e// which are invariant 
under the discrete symmetries of the underlying lattice, but no longer invariant under 
continuous space rotations. Indeed, in the case of the square lattice, the following 
extra term 

2 

h drdrlPd^XP (2.5) 
r=l 
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has to be included in C\j^- Interestingly, in the case of the honeycomb, triangular 
and Kagome lattice, the discrete 60 degrees rotation symmetries do not permit such 
a term - here, it is perfectly legitimate to use the space rotation invariant Lagrangian 
(12 .4p . Note that, for the square lattice, there is an additional contribution with four 
space derivatives 

2 

d r U i d r U l d r U k d r U k . (2.6) 

r=l 

However, as we will show below, terms in that contain four or even more magnon 
fields, are irrelevant for the evaluation of the partition function presented in this work. 

In finite-temperature field theory the partition function is represented as a Eu- 
clidean functional integral 

Tr [exp(-n/T)] = J [dU] exp ( - J d A x £ e// ) . (2.7) 

The integration extends over all magnon field configurations which are periodic in 
the Euclidean time direction U(x, x 4 + (3) = U(x,x^), with f3 = 1/T. The quantity 
C e ff on the right hand side is the Euclidean form of the effective Lagrangian, which 
consists of a string of terms 

C eff = C 2 eff + Ci fJ + 0(p 6 ) , (2.8) 

involving an increasing number of space and time derivatives. 

The virtue of the representation (12. 7p lies in the fact that it can be evaluated per- 
turbatively. To a given order in the low-temperature expansion only a finite number of 
Feynman graphs and only a finite number of effective coupling constants contribute. 
The low-temperature expansion of the partition function is obtained by considering the 
fluctuations of the spontaneous magnetization vector field U = (U 1 , U 2 , U 3 ) around 
the ground state Uq = (0,0,1), i.e., by expanding U 3 in powers of the spin-wave 
fluctuations U a , 

U 3 = Vl - U a U a = 1 - \U a U a - \U a U a U h U b - ... . (2.9) 

Inserting this expansion into formula (12. 7p one then generates the Feynman diagrams 
illustrated in Fig. [TJ The leading contribution in the exponential on the right hand 
side of Eq. (12. 7p is of order p 2 and originates from C> 2 ff. It contains a term quadratic 
in the spin-wave field U a - with the appropriate derivatives and the magnetic field 
displayed in Eq. (l2.2p - and describes free magnons. The corresponding diagram for 
the partition function is the one-loop diagram 4 of Fig. [TJ 

A crucial point underlying the perturbative evaluation of the partition function 
of any system concerns the suppression of loop diagrams in the effective field theory 
framework. In three spatial dimensions - and in the case of the ferromagnet - each 
loop in a Feynman diagram is suppressed by three powers of momentum. In two 



6 




8/ 8g 



Figure 1: Feynman graphs related to the low-temperature expansion of the partition 
function for a two-dimensional ferro magnet up to order p 8 . The numbers attached to 
the vertices refer to the piece of the effective Lagrangian they come from. Vertices 
associated with the leading term C? e ^ are denoted by a dot. Note that ferromagnetic 
loops are suppressed by two powers of momentum in d s =2. 
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dimensions, on the other hand, ferromagnetic loops are only suppressed by two powers 
of momentum: The one- loop diagram 4 is of order p 4 , as it involves C? e ^ (p 2 ) and one 
loop (p 2 ). 

The reason why the loop suppression depends on the spatial dimension d s of the 
system can easily be appreciated as follows: Each loop involves an integral of the type 

dud ds k — — oc p ds , (2.10) 

uj — '-fk 2 

related to ferromagnetic magnons circling in the loop. On dimensional grounds the 
integral is proportional to d s powers of momentum, i.e., each loop in a Feynman 
diagram referring to the two-dimensional ferromagnet is suppressed by p 2 . 

The remainder of the effective Lagrangian in the path integral formula for the par- 
tition function ( 12. 7p . i.e. C i eff + C & eff + . . . , is treated as a perturbation. The Gaussian 



integrals are evaluated in the standard manner (see Ref. |59|, in particular chapter 
3) and one arrives at a set of Feynman rules which differ from the zero-temperature 
rules of the effective Lagrangian method only in one respect: the periodicity condi- 
tion imposed on the magnon fields modifies the propagator. At finite temperature, 
the propagator is given by 

oo 

G(x)= A(f,x 4 + n/3), (2.11) 

n = — oo 

where A(x) is the Euclidean propagator at zero temperature, 

/rlh<rl 2 h p ikx-ik 4 X4 r A^U 

"1 4 , 77 = / 4^ e^ 2 -— . (2.12) 

[2-kY 1 k 2 - ik 4 + fiH J (2tiT 

An explicit representation for the thermal propagator, dimensionally regularized in 
the spatial dimension d s , is 

dg_ OO 2 

G(x) = — Q 2 J2 "IT exp-^-^» B(x n ) , (2.13) 

with 

x n = Xi + nfi . (2.14) 

We restrict ourselves to the infinite volume limit and evaluate the free energy 
density z, defined by 

z = — T lim LT 2 ln[Trexp(--H/T)] . (2.15) 

L— >oo 

Note again that in the case of a quadratic dispersion relation - and in two space 
dimensions - each loop in a Feynman diagram is suppressed by two powers of momen- 
tum. This suppression rule lies at the heart of the organization of the Feynman graphs 
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of the partition function for the two-dimensional ferromagnet depicted in Fig. [TJ Now 
we also understand why terms in that contain four or even more magnon fields 
are irrelevant for the explicit evaluation of the partition function presented in this 
work which goes up to order p 6 : The two- loop diagram 8d with an insertion from 
£>iff, containing four magnon fields, is of order p 8 , as it involves (p 4 ) and two 
loops (p 4 ). 

In the next section we will evaluate the partition function of the two-dimensional 
ideal ferromagnet in full generality up to order p 6 . The evaluation of the partition 
function at order p 8 is much more involved. In particular, the renormalization and 
numerical evaluation of the three-loop graph 8c turns out to be rather elaborate 



- a detailed account of this calculation will be presented elsewhere [60j. Here, we 
rather focus on the general structure of the low-temperature expansion and answer 
the question which contributions originate from noninteracting spin waves and which 
ones are due to the spin-wave interaction - this question has never been addressed so 
far. 



3 Thermodynamics of Two-Dimensional Ideal Fer- 
romagnets 

We now consider those Feynman graphs depicted in Fig. Q] that contribute to the 
partition function up to order p 6 or, equivalently, up to order T 3 . Again, additional 
information on finite-temperature effective Lagrangians and the evaluation of the 
corresponding Feynman diagrams - going beyond the outline given in the previous 
section - can be found in Ref. [9J (see section III and appendix A). 

At leading order (order p 2 ), we have the tree graph 2 involving £} e ^ which merely 
leads to a temperature-independent contribution to the free energy density, 

Z2 = -EfjiH. (3.1) 



The leading temperature-dependent contribution is of order p 4 and stems from the 
one-loop graph 4. It is associated with a (d s + l)-dimensional nonrelativistic free Bose 
gas and amounts to 



47T7 



— T 2 V 



X1 e -fiHnfi 



71=1 



(3.2) 



6a 



At order p 6 the first two-loop graph shows up. This contribution, related to graph 
is proportional to single space derivatives of the propagator at the origin, 



z 6a OC 



d r G(x) 



J x=0 



d r G(x) 



x=0 



(3.3) 
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and thus vanishes because the thermal propagator is invariant under parity, much like 
the Heisenberg Hamiltonian. Remember that the effective Lagrangian - and therefore 
the thermal propagator - inherits all the symmetries of the underlying Heisenberg 
model. 

At the same order p 6 , the next-to-leading order Lagrangian comes into play 
The one-loop graph 6b, which involves a two-magnon vertex, corresponds to 



Z6b — —- 



2h 



A 2 G(x) 



2L 



J x=0 



r=l 



J x=0 



(3.4) 



and yields the temperature-dependent contribution 

T 4/ 3 + 3/4 



°6b 



47rS7 3 



^ e -fiHnf3 



71=1 



7V 



(3.5) 



Collecting terms, the free energy density of the two-dimensional ideal ferromagnet 
up to order p 6 cx T 3 becomes 

^ 47T7 n 2 47rS7 3 n 3 ^ ^ ^ ^ 

' 71=1 ' 71 = 1 

The contributions of order T 2 and T 3 arise from one-loop graphs and are related 
to the free energy density of noninteracting spin waves. The former contribution is 
exclusively determined by the leading-order effective constants £ and F (7 = _F 2 /E), 
i.e., it is the same for any of the four types of lattices - square, honeycomb, triangular 
and Kagome - considered here. The term of order T 3 , on the other hand, is not 
universal since it involves the next-to-leading order effective constants Z3 and 1 4 from 
l~-eff- In the case of the honeycomb, triangular and Kagome lattice, the coefficient of 
order T 3 exclusively contains the contribution from Z 3 - the term (12.51) in the effective 
Lagrangian involving Z4, which accounts for the lattice anisotropies, is excluded due 
to the discrete 60 degrees rotation symmetries of these lattices. Remarkably, the 
spin-wave interaction does not yet manifest itself at this order in the low-temperature 
expansion of the free energy density. The only potential candidate, the two-loop 
diagram 6a of order T 3 , turns out to be zero due to parity. 

The ratio fiH(3 = [iH/T in the above series can take any value, as long as the 
temperature and the magnetic field themselves are small compared to the intrinsic 
scale of the underlying theory, which in the present case of the two-dimensional ferro- 
magnet is given by the exchange integral J of the Heisenberg model. In the following 
we will be interested in the limit T 3> [iH which we implement by holding T fixed 
and sending the magnetic field to zero. Since we keep the fixed temperature small 
compared to the scale J, we never leave the domain of validity of the low-temperature 
expansion. 
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In order to discuss the effect of a weak magnetic field we thus expand the result 
(13.61) in the dimensionless parameter 



a = ^H(3 = . (3.7) 



T 

Retaining all terms up to quadratic in a, we obtain 

z = -S/iiJ-^-T 2 |c(2) + (7lna-a- ^- + 0(a 3 )} 
4Z 3 + 3L , f >/0 v 1 2l 3a 2 



47rS7 3 



T 3 {C(3) - C(2) a - V In a + ^- + 0(a 3 ) } + (9(p 8 ) . (3.8) 



In the absence of an external magnetic field, the sums in the series (I3.6P become 
temperature independent and reduce to Riemann zeta functions, 

' C(2)T 2 -^±^C(3)T 3 + 0(/). (3.9) 



47T7 47rS7 

Since we are dealing with a two-dimensional system, we have to be careful by tak- 
ing the limit H — > due to the Mermin- Wagner theorem. A thorough discussion, 
confirming the validity of the above series, will be given in Sec. HI 

To corroborate the structure of the low-temperature series, let us consider an 
independent derivation, based on the evaluation of the two-point function and the 
subsequent extraction of the dispersion relation. The free energy density of a two- 
dimensional gas of noninteracting bosons is then obtained from the dispersion relation 
through 

z = z + — [d 2 k In [l - e -^ T ] , (3.10) 



(2vr) 2 j 

where zq is the free energy density of the vacuum. The leading term in the dispersion 
relation, 

u(k)= 1 k 2 + ^H + 0(\k\ A ) ) 7= If, ( 3 - U ) 

yields the dominant one-loop contribution zj oc T 2 in the free energy density of the 
two-dimensional ferromagnet. 

Subleading terms in the dispersion relation are obtained by evaluating the two- 
point function to higher orders. The relevant graphs are shown in Fig. [2J Depicted are 
all contributions up to order p 6 . Instead of listing individual results for the two-point 
function, we give the final expression for the dispersion relation originating from these 
graphs: 

u (k) = I jE/iif + F 2 k 2 - (2/ 3 + lh)k 4 + C(fc 6 )} . (3.12) 

Note that the one-loop graph 6a does not contribute to the dispersion relation. There 
is thus only one additional diagram contributing at the order we are considering: 
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4 



6 a 



66 



Figure 2: Feynman graphs occurring in the low-energy expansion of the two-point 
function for a two-dimensional ferromagnet up to order p 6 . 



Graph 6b which leads to a higher-order term involving the effective constants I3 and Z4. 
Again, in the case of the honeycomb, triangular, and Kagome lattices, the contribution 
proportional to the low-energy constant I4 is absent. Inserting the above expression 
into the free Bose gas formula (13.101) one readily confirms the low-temperature series 
for the free energy density. 

Let us also consider the low-temperature series for the energy density u, for the en- 
tropy density s, and for the heat capacity cy of the two-dimensional ideal ferromagnet. 
They are readily worked out from the thermodynamic relations 

dP „ ^ du ds 

s — -z=; , u = Ts — P, Cy = — - = T — — . (3.13) 
dT ' dT dT K J 

Because the system is homogeneous, the pressure can be obtained from the temperature- 
dependent part of the free energy density, 

P = z Q -z, (3.14) 

such that the other thermodynamic quantities amount to 



u 



1 I 00 „-cm 00 

-T 2 aV — + E 



47T7 ^— ' n ^— ' n 

' n=l n=l 



2 



00 _„„ 00 



' I n=l n=l 



00 -cm 



' I n=l n=l 



' I n=l n=l 



,— cm 



00 -crn 



cy = J-T a 2 Ve— + 2aVE_ + 2 V 



2 



V n=l n=l n=l 

I n=l n=l n=l 
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For a weak magnetic field H, the series may be expanded in the parameter a = fj,H/T, 



u 



Li* { f(2) _, + £ + «,(„.)} 
^±pir»{2 C (3)-c(2). + ^ + oM} + o^ 



— T \2((2) + a In a - 2a + Ola 3 ) } 

+ 7,2 ( 3C(3) - 2C(2)ff ~y ln<r+ ^r + °' ff4 )} + °<f 6 



2 



i-T{2C(2)-a + (9(a 3 )} 

T 2 {6C(3) - 2C(2)a + ^ + + 0(p 6 ) , (3.16) 



where we have retained terms up to quadratic in the magnetic field. Formally, as it 
was the case for the free energy density, the limit H — > poses no problems. Note 
again that all contributions in the above series for u, s and cy originate from one-loop 
graphs - the spin-wave interaction does not yet manifest itself at this order of the 
low-temperature expansion. 

We may now compare our results, derived within the effective field theory frame- 
work, with the literature. The few explicit results available all refer to the limit H — >0 
and to the square lattice 10l.ll3l4l7l|. For the free energy density of the two-dimensional 
ideal ferromagnet these authors obtain 

C(2)T 2 - l— 9 C(3)T 3 + 0(T 4 ). (3.17) 



Art J So 2 * v ' 32tt J 2 S 2 a 2 

The first term coincides with our series (I3.9P provided that we express the effective 
low-energy constant 7 in terms of microscopic constants as 

7 = JSa 2 . (3.18) 

By matching the coefficients of the second term, the combination of the low-energy 
constants l 3 and 1 4 is identified as 

o , J S GL /„\ 

h + ^h = ^^. (3.19) 

The microscopic calculation is thus consistent with the effective calculation for the 
square lattice in the limit H— >0. 



However, the author of Ref. [13j, which is considered as a standard reference on 
the low-temperature properties of two-dimensional ferromagnets, was rather cautious 
about the correctness or validity of his result: Regarding the low-temperature series 
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for the free energy density he comments that it is possible that his series gives the 
correct low-temperature expansion. 



The reason for his caution may be readily identified. While spin-wave theory 
works well for three-dimensional systems, in two or one space dimensions the spin- 
wave expansion is plagued with divergences. In order to cope with low-dimensional 
systems, many approximations were invented. One very popular method is modified 



spin-wave theory, advocated first for two-dimensional Heisenberg ferromagnets [13 



and then transferred to two-dimensional antiferromagnets [6lH63| . The essential idea 
is to impose an ad hoc condition on the chemical potential. However, to the best 
of our knowledge, the justification of such an ad hoc condition was never rigorously 
examined, neither for the ferromagnet nor for the antiferromagnet. 

Within our effective field theory framework, we will put the above low-temperature 
series for a two-dimensional system on a firm basis - on the same footing as the low- 
temperature series for ferro- and antiferromagnets in three space dimensions. Indeed, 
as we will show in the next section, the Mermin- Wagner theorem is perfectly con- 
sistent with the low-temperature series derived in the present work. We also like to 
emphasize that the effective field theory approach we have pursued does not resort 
to any approximations or ad hoc conditions as e.g. in the case of modified spin-wave 
theory. Moreover, our series go beyond the results of the literature, as they explicitly 
include a weak external magnetic field and are valid not only for the square lattice, but 
also for the honeycomb, the triangular and the Kagome lattice with a spontaneously 
broken spin symmetry 0(3) — > 0(2). 



4 Range of Validity of the Low- Temperature Series 

Whereas in three space dimensions the limit can readily be taken, one has 



to be careful in two (or one) dimensions due to the Mermin- Wagner theorem [64J, 
which states that no spontaneous symmetry breaking at any finite temperature in the 
0(3)-invariant two-dimensional Heisenberg model can occur. Accordingly, no mass- 
less magnons in the low-energy spectrum at any finite temperature will be present. 
In the context of ferromagnetic magnons this means that the low-energy spectrum 
exhibits a nonperturbatively generated energy gap and that the correlation length of 
the magnons no longer is infinite. Still, the correlation length is exponentially large, 



the argument of the exponential being proportional to the inverse temperature [65 



i / T 



£ np = Ct.aS f \/j^exp . (4.1 



T 



Here a is the spacing between two neighboring sites on the square lattice and the 
quantity ~ 0.05 is a dimensionless constant. 

Strictly speaking, it is therefore not legitimate to switch off the external magnetic 
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field H completely, because our effective calculation does not take into account the 
nonperturbative effect. However, the corrections due to the nonperturbatively gen- 
erated energy gap are so tiny, that they cannot manifest themselves in the power 
series derived in this work. In what follows, we will estimate the order of magnitude 
of these corrections and thus verify this claim. While the above explicit expression 
for the correlation length refers to the square lattice, note that the conclusions to be 
presented in this section also apply to the honeycomb, the triangular and the Kagome 
lattice. 

Our low-temperature series are valid as long as the correlation length £ of the 
Goldstone bosons is much smaller than the nonperturbatively generated correlation 
length £ np . In order to define the correlation length £ for ferromagnetic magnons in 
a natural way, let us consider the dispersion relation. In the presence of a magnetic 
field it takes the form 

u(k)= 1 k 2 + f iH + 0(\k\ A ), 7 =S, (4-2) 



S 



and we may define the correlation length as 



This quantity has dimension of length and tends to infinity for zero magnetic field, as 
it should. It is the analog of the corresponding relation for antiferromagnetic magnons, 
which obey a linear (relativisticjdispersion law. In that case the correlation length is 
given by the inverse mass M n [36], 

where S s and H s are the staggered magnetization at zero temperature and the stag- 
gered field, respectively. 

The low-temperature series derived in the previous section are certainly valid if 
£ np is - let us say - a thousand times larger than £, i.e. 



1 g 

1000 i np 






r 2ttJS 2 i 


Lh exp 


I T \ 



(4.5) 



Note that we have used eq. fl 3 . X 8 j) in order to express the effective constant F in terms 
of the exchange integral J of the underlying theory by 



F = VZJSa . (4.6) 

Now the exchange integral defines a scale in the underlying theory and for the effective 
expansion to be consistent, the temperature has to be small with respect to this scale. 
Assuming that 

7=i5o< < 4 - 7 ' 
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relation (14. 5p then yields the ratio 



^1(T 125 (£ = I). (I. 



We thus see that, in principle, we cannot completely switch off the magnetic field - 
rather, we start running into trouble as soon as the ratio fiH/T is of the order of the 
above value. However, the error introduced into the low-temperature series considered 
in this work is indeed extremely small. Hence we confirm that the corrections due 
to the nonperturbatively generated energy gap are so tiny that they cannot manifest 
themselves in the above low-temperature expansions for the thermodynamic quantities 
- the subtleties raised by the Mermin- Wagner theorem in d s =2 are not relevant for 
our calculation. 

The effective calculation performed in this work is restricted to the regime £ <C £ np . 
This does not mean that the regime £ 3> £ n p is beyond the reach of the effective field 
theory. Rather, one has to resort to a different type of perturbative expansion. A 
similar situation arises when one considers finite-size effects. In particular, when the 
Goldstone boson mass is small compared to the inverse size of the box, a different ef- 
fective expansion scheme, the so-called e-expansion, applies. Indeed, various problems 



within this framework have been investigated in detail [661-171 



We close this section with a conceptual remark. Our effective analysis refers to 
the two-dimensional ideal ferro magnet, i.e., to a two-dimensional system which is gov- 
erned by the isotropic exchange interaction and the interaction with a weak external 
magnetic field. This represents the system which was analyzed before within a micro- 



scopic framework by Takahashi and other authors |10Hl9| . For this idealized situation 
we have rigorously shown that taking the limit H — > in the low-temperature series 
derived in this work is consistent with the Mermin- Wagner theorem. 

In a more realistic approach to ferromagnetic films, however, one has to also con- 
sider the magnetic anisotropy and dipolar interactions. Although they are much 
weaker than the exchange interaction, these effects may play a decisive role (see 



e.g. 72|). In particular, taking into account these effects, the Mermin- Wagner theorem 
is evaded since two of its basic assumptions are no longer fulfilled: the Hamiltonian 
is no longer isotropic and the dipolar interaction is no longer short-ranged. 

So the question regarding the implications of the Mermin- Wagner theorem on the 
low-temperature properties of a two-dimensional ideal ferromagnet is rather academic. 



Still, unlike the various authors before |1CH19|. here we have put the low-temperature 



series on a firm basis in this idealized framework. 
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5 Ideal Ferromagnets in Two and Three Spatial 
Dimensions — A Comparison 



It is very instructive to compare the thermodynamic properties of two-dimensional 
ferromagnets with those of three-dimensional ideal ferromagnets within the effective 
field theory perspective, adopting thereby a unified point of view based on the sym- 
metries of the system. 

As we pointed out in Sec. [2J the suppression of loops in the perturbative expansion 
of the partition function depends on the spatial dimension. For two-dimensional ferro- 
magnetic systems, loops are suppressed by two powers of momentum, in three spatial 
dimensions, on the other hand, loops are suppressed by three powers of momentum. 
Accordingly, the organization of Feynman diagrams related to the three-dimensional 
ferromagnet, depicted in Fig. [3] (for details see Ref. jl]), is quite different from the one 
referring to the two-dimensional ferromagnet, Fig. [TJ Still, as we now discuss, there 
are also various similarities between d s =2 and d s =3. 

In either case the leading temperature-dependent contribution to the free energy 
density is related to a one-loop graph. For the three-dimensional ideal ferromagnet it 
is of order p 5 oc T 5 / 2 , for the two-dimensional ideal ferromagnet we have p 4 oc T 2 . The 
next-to-leading contribution again stems from a one-loop graph, but this time with 
an insertion from £ 4 yj. For the three-dimensional ferromagnet, this term - diagram 
7 of Fig. 0- is of order p 7 oc T 7//2 . For the two-dimensional ferromagnet, diagram 6b 
of Fig. [T] leads to a term of order p 6 oc T 3 . 

Beyond one-loop order the spin-wave interaction comes into play. However, the 
corresponding two-loop graph (graph 6a for d s =2, graph 8 for d s =3) which would 
represent the first candidate for the spin-wave interaction, turns out to be zero due 
to parity. Regarding the three-dimensional ideal ferromagnet there is thus no contri- 
bution of order p 8 oc T 4 in the series for the free energy density. Likewise, for the 
two-dimensional ideal ferromagnet this means that the contribution of order p 6 oc 
in the free energy density is exclusively related to noninteracting spin waves. 

So at which order does the spin-wave interaction manifest itself in either case? 
In three dimensions it shows up at order p 10 oc T 5 due to the two-loop diagrams 
which involve insertions from the next-to-leading Lagrangian £ 4 yj, i.e. graphs 10a 
and 10b of Fig. El In the case of the two-dimensional ideal ferromagnet the spin-wave 
interaction sets in at order p 8 oc T 4 through five different diagrams - two-loop as well 
as three-loop diagrams (see Fig. [1]). Note that for the three-dimensional ferromagnet 
these three- loop diagrams are of order p 11 oc T 11 / 2 , i.e. of higher order than the two- 
loop diagrams 10a and 10b with insertions from £ 4 yj. In two spatial dimensions, on 
the other hand, they are of the same order as the two-loop diagrams, all five graphs 
contributing at the same order p 8 oc T 4 . The explicit evaluation of these contributions 



is quite involved and will be presented elsewhere 60] - here we rather want to draw 



17 




lid lie 



Figure 3: Feynman graphs related to the low-temperature expansion of the partition 
function for a three-dimensional ferromagnet up to order p 11 . The numbers 
attached to the vertices refer to the piece of the effective Lagrangian they come from. 
Vertices associated with the leading term C? e ^ are denoted by a dot. Note that 
ferromagnetic loops are suppressed by three powers of momentum in d s —3. 
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our attention to the general structure of the low-temperature series. 



Still, we have to mention that for the two-dimensional ideal ferromagnet, unlike 
for the ideal ferromagnet in three spatial dimensions, logarithmic contributions in 
the low-temperature series will show up. This is related to the structure of the ul- 
traviolet divergences arising in higher-order loop-diagrams, which in the case of the 
two-dimensional ferromagnet require a logarithmic renormalization of next-to-leading 



order effective coupling constants - again, details will be provided in Ref. [60 



Summarizing the above results, the low-temperature expansion for the free energy 
density of the ideal ferromagnet in two and three spatial dimensions - in the absence 
of an external magnetic field - exhibits the following general structure, 



z da =2 = -VoT'-v^ + OiTHnT,^), 

z ds=3 = —h T^ — h\T^ — h 2 T% — h 3 T 5 + C(T 11//2 ) , (5.1) 
where we have highlighted all terms which are related to the spin-wave interaction. 

Note that in the series for the two-dimensional ideal ferromagnet no half-integer 
powers of the temperature occur. The first contribution is of order T 2 and any other 
corrections necessarily involve integer powers of the temperature. This is because each 
additional loop yields an additional power of T. Likewise, higher-order vertices with 
insertions from the effective Lagrangian, 

C eff = C 2 eff + C\ fJ + C% f + 0(p 8 ) , (5.2) 

also increase the temperature power in steps of p 2 ocT. 

Now in three dimensions, the first contribution in the free energy density is of 
order T 5 / 2 . Also here, insertions of higher-order contributions from the effective La- 
grangian yield additional integer powers of the temperature. Successive insertions of 
higher-order vertices in one-loop graphs, e.g., lead to the pattern T 7 / 2 , T 9//2 , T 11//2 , . . . , 
describing the effect of noninteracting spin waves. However, since loops in three spa- 
tial dimensions are suppressed by three powers of momentum, or equivalently, lead to 
additional powers of T 3//2 , we will also have integer powers of the temperature in the 
above series. In fact, any such integer power in the series for the three-dimensional 
ferromagnet necessarily must have its origin in the spin-wave interaction. 



6 Conclusions 



The present study was devoted to the thermodynamic properties of two-dimensional 
ideal ferromagnets at low temperatures. Previous articles, based on modified spin- 
wave theory or Schwinger Boson mean field theory, have also discussed the low- 
temperature behavior of two-dimensional ferromagnets, but there the magnons were 
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considered as ideal Bose particles - the problem of the spin-wave interaction was 
neglected and it thus remained unclear whether the low-temperature series given in 
these articles are complete or will receive additional corrections due to the interaction. 
Furthermore, it has never been discussed in a systematic manner how a weak external 
magnetic field manifest itself in the low-temperature series for the thermodynamic 
quantities or how these series look like for underlying geometries other than a square 
lattice. 

Within the effective Lagrangian framework, we have addressed all these questions 
in detail. We have derived the low-temperature expansion of the partition function 
up to two-loop order - i.e. order T 3 - for two-dimensional ideal ferromagnets on 
the square, honeycomb, triangular and Kagome lattice, where the 0(3) spin rotation 
symmetry is spontaneously broken to 0(2) by the ground state. Remarkably, the spin- 
wave interaction does not yet manifest itself at order T 3 in the free energy density - 
it will only enter at order T 4 InT. Analogously, in the case of the three-dimensional 
ideal ferro magnet, the spin- wave interaction does not yet manifest itself at order T 4 
- rather, as Dyson showed, it enters at order T 5 in the free energy density. In both 
cases the spin-wave interaction is thus very weak. 

While the validity of the low-temperature series derived within the framework of 
modified spin-wave theory - which resorts to an ad hoc assumption - has never been 
critically examined, here we have put these series on safe grounds. Indeed, as discussed 
in detail, the Mermin- Wagner theorem is perfectly consistent with our results and the 
low-temperature series are valid as they stand. 

In conclusion, the effective field theory method is a very powerful tool to analyze 
the general structure of the low-temperature expansion of the partition function for 
systems exhibiting collective magnetic behavior. Not only have we conclusively dis- 
cussed the effect of the spin-wave interaction and a weak magnetic field in a systematic 
manner, but also have we put our results on a firm basis. In a more realistic approach 
to ferromagnetic films, one should also include the magnetic anisotropy and dipolar 
interactions. Here, much like Dyson, Takahashi and other authors, we have considered 
the ideal ferromagnet and rigorously derived the low-temperature properties for this 
"clean" system. 
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